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RESUMEN 

Actualmente esta en construccion ocho bloques estructurales, donde 
funcionaran los Centros de Investigation, Postgrado y Administrative de la 
Universidad de Fuerzas Armadas ESPE. Los bloques: 1, 2, 3, 6, 7 y 8 tienen 
aisladores sfsmicos de triple pendulo de friccion y los bloques: 4 y 5 tienen 
disipadores de energia. 

En este articulo se realiza el analisis sfsmico espacial del bloque estructural 
8, considerando tres grados de libertad por el Metodo Espectral; y luego se realiza 
un analisis sfsmico piano de dos porticos, con el proposito de comparar resultados. 

Se utiliza el modelo de tres fases de Me. Vitty y Constantinou (2015) para 
definir la rigidez efectiva del elemento aislador, para el sismo de diseno DBE y para 
el sismo maximo considerado MCE. 

Palabras claves. - Aislador de triple pendulo de friccion, metodo modal espectral, 
maximo sismo considerado. 



ABSTRACT 

Currently, eight buildings are under construction in the Armed Forces 
University - ESPE. There will function the Graduate Investigation Center and 
Administrative Department. The buildings 1 , 2, 3, 6, 7 and 8 are being constructed 
with seismic isolators (triple frictional pendulum), while the buildings 4 and 5 are 
being constructed with energy dissipaters. 

This article shows seismic analysis by different techniques in the building 8 
which are composed for frames 44 and 48. Moreover, it makes a comparison the 
dynamic properties between a 3D analysis for the building 8 and 2D analysis for the 
frames 44 and 48. It considers three freedom degrees per floor, two components of 
horizontal displacement and one rotation per floor. 
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These analyses are conducting with Spectral Modal Method which are 
composed by three phases; to define the effective stiffness of the isolate element, 
develop design basis earthquake and determine maximum credible earthquake. 

Key words: seismic isolators (triple frictional pendulum), Spectral Modal 
Method and maximum credible earthquake. 

1. INTRODUCTION 

El Bloque Estructural 8 de la UFA-ESPE, esta destinado a la parte 
Administrativa de los Post Grados, tiene dos pisos y los aisladores se van a colocar 
sobre los pedestales que sobresalen de la cimentacion. 

En la figura 1 se presentan varias fotograffas tomadas las de arriba el 12 de 
diciembre y las de abajo el 28 de diciembre de 2015, en ellas se aprecia el replantillo 
de una cadena; la fundicion de un plinto (las fotos superiores) y el armado de las 
cadenas cuya armadura longitudinal atraviesa los pedestales sobre los que van los 
aisladores sfsmicos (las dos inferiores). 





Figura 1 Panoramica del Bloque Estructural 8. 
(Fotograffas tomadas en diciembre de 2015). 
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Al colocar los aisladores sobre los pedestales de los plintos se debe 
construir una losa de aislacion, sobre ellos. De tal manera que la estructura tendra 
tres losas en total, a saber: la de aislacion (cuyas vigas son las mismas del primer 
piso), la del primer piso y la del segundo piso, las mismas que se muestran en la 
figura 2, con sus respectivos niveles. 

En la figura 2 se observa la parte circular de la losa, a la que hace referenda 
la fotograffa izquierda de la figura 1 y la parte central donde cambia de direccion la 
estructura. Las losas de los niveles 4.94 y 9.88 tienen voladizos los mismos que no 
se dibujan en los modelos, pero si son considerados en el analisis sfsmico. 

Se observa que la estructura es irregular en planta, por lo que el analisis 
sfsmico se realiza de dos maneras; la primera espacial y la segunda plana. 




Nivel +0.04 



Nivel +4.94 




Nivel +9.88 



Figura 2 Losas del Bloque Estructural 8 de la UFA-ESPE. 



2. ANALISIS SISMICO ESPACIAL 

En la figura 3 se presenta el modelo de calculo. Basicamente no se ha 
considerado el cambio de direccion de la planta para que los porticos longitudinales 
sean rectos. 




Figura 3 Esquema del Modelo Adoptado 
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Figura 4 Losa Nivel 0.04, superior izquierda; Nivel 4.94, superior Derecha; Nivel 

9.88, Inferior 
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2.1 Centro de Gravedad de cada planta 

En la tabla 1 se resume el calculo del centra de gravedad de cada una de 
las plantas de acuerdo al modelo de calculo adoptado; se ha considerado como 
origen de coordenadas, la intersection del eje 44 con el eje Y (Parte inferior 
izquierda), los Centros de Gravedad encontrados se indican en las figuras 5. 



Tabla 1 Centro de Gravedad de las losas del Modelo Numerico adoptado 



Nivel 




Ycg 


0.04 


6.25 


16.99 


4.94 


6.29 


17.00 


9.88 


6.50 


16.88 




Figuras 5 Centros de gravedad de las losas del Modelo Numerico adoptado 



Se considera que las cargas se encuentran distribuidas en forma simetrica. 
En este caso el centra de gravedad coincide con el centra de masas. 



2.2 Cargas que Megan a los aisladores 



El Bloque Estructural 8 tiene 9 aisladores, que han sido identificados como 
se indica en la grafica superior izquierda de la figura 6. En las restantes graficas se 
indica el area cooperante que gravita en cada uno de los aisladores; en base a esta 
area se determino el peso que llega a cada aislador; para el efecto se trabajo con el 
siguiente porcentaje de carga viva: 

CT = D + 0.25 L O) 



Tabla 2 Cargas de las losas del 


Modelo Numerico adoptado 


Carga Muerta 


D = 0.8 7„ 2 


Losa de aislacion y primer piso 


O = 0.6 r/ m2 


Losa de cubierta (segundo piso) 


Carga Viva 


L = °- 25 r / m 2 


Todas las losas 
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En las figuras 6 se observa que se han definido tres ejes longitudinales para 
determinar la posicion de los aisladores y son los denominados Y, Y’, Z. 






Figuras 6 Areas cooperantes para el calculo de la carga que Mega a los 

aisladores. 
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A manera de ejemplo se indica en las tablas 3 y 4 el calculo de las cargas 
que Megan a los aisladores 44 Y, 44 Z. 



Tabla 3 Calculo de la Carga que llega al aislador 44Y 



Aislador 44Y 


Nivel 


Area Cooperante 


Carga Vertical 


W 


(m 2 ) 


(T/m 2 ) 


(T) 


+ 0.04 


57.97 


0.86 


50.00 


+ 4.94 


82.55 


0.86 


71.20 


+ 9.88 


82.55 


0.66 


54.69 




Total 


175.88 



Tabla 4 Calculo de la Carga que llega al aislador 44Z 



Aislador 44Y 


Nivel 


Area Cooperante 
(m 2 ) 


Carga Vertical 
(T/m 2 ) 


W 

(T) 


+ 0.04 


30.47 


0.86 


26.28 


+ 4.94 


52.58 


0.86 


45.35 


+ 9.88 


52.58 


0.66 


34.83 




Total 


175.88 



2.3 Rigidez efectiva de los aisladores 

En la figura 7 se presenta el aislador FPT8833/1 2-1 2/8-6, que se colocara 
en el Bloque Estructural 8. En base a la geometna de este aislador y considerando 
los coeficientes de race del aislador interior: n 2 = n 3 = 0.01 ; y para las placas 
exteriores: ^ = ju 4 = 0.07, y las cargas verticales que Megan al aislador se determino 
la rigidez efectiva para los sismos DBE y MCE utilizando el modelo de 3 fases 
propuesto por Me Vitty y Constantinou (201 5). Los coeficientes de amortiguamiento 
equivalente para el sismo DBE es % eq = 0.2028 y para el sismo MCE % eq =0.1962. 




Figura 7 Aislador FPT8833/1 2-1 2/8-6, utilizado en el Proyecto. 



A manera de ejemplo se indica en las tablas 5, 6 y 7 el calculo de rigideces 
efectivas para los aisladores de los porticos Y, Z, inclinado izquierdo e inclinado 
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derecho, considerando a estos dos ultimos que parten del eje 48 y llegan al eje Y 
como del eje Z. 



Tabla 5 Calculo de la rigidez efectiva de los aisladores Portico Y Bloque 8 



Tabla resumen del Portico Y Bloque 8 


Aislador 


w 


Fmax 2 


q2 


q3 


Kef2 


Kef3 


(T) 


(T) 


(m) 


(m) 


(T/m) 


(T/m) 


44Y 


175.88 


33.299 






62.056 


58.801 


47Y 


297.22 


56.272 


0.537 


0.566 


104.868 


99.368 


48Y 


161.82 


30.637 






57.095 


54.101 



Tabla 6 Calculo de la rigidez efectiva de los aisladores Portico Z Bloque 8 



Tabla resumen del Portico Z Bloque 8 


Aislador 


W 


Fmax 2 


q2 


q3 


Kef2 


Kef3 


(T) 


(T) 


(m) 


(m) 


(T/m) 


(T/m) 


44Z 


106.46 


20.156 






37.562 


35.592 


45Z 


147.63 


27.951 






52.088 


49.357 


46Z 


79.96 


15.139 


0.537 


0.566 


28.212 


26.733 


47Z 


153.21 


29.007 






54.057 


51.222 


48Z 


161.82 


30.637 






57.095 


54.101 



Tabla 7 Calculo de la rigidez efectiva de los aisladores Portico Inclinado 
Derecho / Izquierdo Bloque 8 



Tabla resumen del Portico Inclinado Derecho / Izquierdo Bloque 8 


Aislador 


w 

(T) 


Fmax 2 
(T) 


q2 

(m) 


q3 

(m) 


Kef2 

(T/m) 


Kef3 

(T/m) 


49Y' 


83.00 


15.715 


0.537 


0.566 


29.285 


27.749 


48Y / 48Z 


161.88 


30.649 


57.116 


54.120 



2.4 Distancia del centro de masas a los porticos inclinados 

Para determinar la matriz de compatibilidad de deformaciones A, que 
relaciona las coordenadas de piso con las coordenadas laterales de los porticos, se 
necesita encontrar la distancia del centro de masas a los porticos. Cuando se tiene 
estructuras por porticos ortogonales el calculo es directo, pero cuando se tienen 
porticos inclinados se debe proceder como se indica a continuation para la losa del 
sistema de aislacion. 
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Figura 8 Calculo de la distancia del Centro de Masa a los porticos inclinados. 



Con relacion al Centro de Masas (CG) se ha determinado las coordenadas 
de los puntos donde se encuentran los aisladores, de tal manera que se tienen 2 
puntos con sus respectivas coordenadas, lo que permite calcular las ecuaciones de 
las rectas denominadas izquierda (LI) y derecha (L2). 



y-0.96x- 19.91 = 0 


(LI) 


(2) 


y + 0.96x — 19.91 = 0 


(L2) 


(3) 



La distancia de un punto de coordenadas X 0 ,Y 0 -, a una recta de la forma 
Ay + Bx + C = 0, viene dada por la siguiente ecuacion. 



\Ay 0 + Bxq + C| 
y[A 2 + B 2 



( 4 ) 



Luego las distancias del Centro de Masas a las rectas LI vale: 



d 1 = 



\Ay 0 + Bx 0 + C | 
VA 2 + B 2 



|-19.91| 
Vl 2 + (-0.96) 2 



14.363 
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La distancia a la otra recta vale lo mismo, de igual manera se procede con 
los restantes pisos. Las distancias del centra de masas a los porticos se indica en 
la tabla 8, con el signo que tienen: positivo si la orientacion positiva del portico rota 
en sentido anti horario con respecto al centra de masas. (Aguiar 2012, 2014). 



Tabla 8 Distancias del Centro de Masas a los porticos 



Portico 


Distancia d (m) 


Nivel +0.04 


Nivel +4.94 


Nivel +9.88 


Y 


-6.25 


-6.289 


-6.53 


Z 


6.25 


6.211 


5.97 


44 


16.994 


17.883 


17.723 


47 


1.556 


0.667 


0.827 


48 


-13.906 


-13.017 


-13.177 


Incl. Derecha 


14.36 


13.746 


14.028 


Incl. Izquierda 


-14.36 


-13.692 


-13.645 



2.5 Geometria de los porticos 

Las columnas del Bloque 8 son cuadradas, tubulares de 550/550/20 mm, 
las vigas son tipo “I” de 550/300/25/12.5 mm, estas son las vigas principales y para 
el analisis sismico se considera que solo estas vigas resisten la action sismica, 
modelo conservador, debido a que existen vigas secundarias, terciarias y 
cuaternarias. En la figura 9 se indica la geometria de las vigas principales y 
columnas de acero. 



S«pwi 
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o| 
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00424 
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Figura 9 Secciones transversales de los elementos viga y columna del Bloque 8. 

En la figura 1 0 se muestra la geometria de todos los porticos que conforman 
el modelo numerico de calculo adoptado. En la parte central se indica la geometria 
en planta y en la periferia se presentan los porticos de este modelo. 
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tt ODBJOd 

Figura 10 Porticos del Bloque 8 que contempla el Modelo adoptado 



2.6 Matriz de compatibilidad de deformaciones A 

La matriz de rigidez en coordenadas de piso K E se encuentra con la 
siguiente ecuacion. (Aguiar 2012. 2014) 

K E = 'E,A® t K® A® ( 5 ) 



Donde K E es la matriz de rigidez en coordenadas de piso; K L es la matriz 
de rigidez lateral; 4 (i) es la matriz de compatibilidad. que tiene la siguiente forma 
para los porticos del bloque estructural 8. La sumatoria se extiende a todos los 
porticos. 





"cos 6 


0 


0 


send 


0 


0 di 


0 


0 


A& = 


0 


cosd 


0 


0 


send 


0 o 


d 2 


0 




. 0 


0 


cosd 


0 


0 


send 0 


0 


d-3- 



Donde 9 es el angulo que forma la orientacion positiva del portico con el eje 
de las X. Para los porticos paralelos al eje X (Porticos 44, 47 y 48) el valor de 9 = 0; 
y para los porticos perpendiculares al eje X (Porticos Y, Z) 9 = 90. Para el portico 
inclinado izquierdo 9 = 45; para el portico inclinado derecho 9 = 135. A manera de 
ejemplo se indican las matrices de compatibilidad de los dos porticos inclinados. 
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AilnclD ) _ 



-0.707 

0 

. 0 



0 



-0.707 

0 



0.707 0 

0 0.707 

-0.707 0 



0 0 14.360 

0 0 0 

0.707 0 0 



0 0 ' 
13.746 0 

0 14.028. 



Iz ) _ 



0.707 

0 

0 



0 0.707 0 0 0-14.360 

0.707 0 0.707 0 0 0 

0 0.707 0 0.707 0 0 



0 0 
-13.746 0 

0 -14.028. 



2.7 Matriz de Masas 

En la tabla 9 se presenta el area, peso y masa total de cada uno de los pisos 
de acuerdo al modelo de diseno adoptado. 

Tabla 9 Determinacion de las masas de cada planta 



Nivel 


Area 


Peso 


Masa 


(m 2 ) 


(T) 


(T*s 2 /m) 


0+0.04 


423.75 


365.48 


37.26 


0+4.94 


648.06 


558.95 


56.98 


0+9.88 


668.15 


442.65 


45.12 



Para el calculo del momento de inercia de la masa J se divide las plantas de 
las losas en dos figuras, una rectangular y otra triangular y se halla el momento de 
inercia con respecto al centra de masas, aplicando el teorema de los ejes paralelos. 

Jcm = ZJcMi + df ( 6 ) 

Donde ] CMi es el momento de inercia de la masa con respecto al centra de 
masas de la losa rectangular o triangular; es la masa de la losa rectangular o 
triangular y d; es la distancia desde el centra de masas de la losa rectangular o 
triangular al centra de masas de toda la planta. Los valores del momento de inercia 
encontrados son los siguientes: ] 1 = 4053.09 T m s 2 ; J 2 = 7717.68 T ms 2 -, J 3 = 
5940.41 T ms 2 



2.8 Espectros de analisis 

Para el Modelo Numerico de calculo, adoptado para el analisis sfsmico 
espacial y piano, se considera que la estructura es regular en planta y elevacion por 
lo que 0 P = 1; 0 e = 1. 

Por otra parte, se trabaja con un factor de reduccion de las fuerzas sismicas 
R = 1, de tal manera que la estructura va a trabajar en el rango elastico, son los 
aisladores los que van a disipar la energfa sfsmica. 
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El espectro del sismo de diseno DBE es el que reporta el NEC-15 para 
Sangolqiu en un perfil de suelo C; y el espectro del sismo maxi mo considerado MCE 
se halla multiplicando el sismo de diseno por 1.4. Constantinou etal. (2016). 

En la figura 11 se presentan estos espectros, para el sismo DBE a la 
izquierda y para el MCE a la derecha. La Ifnea azul corresponde al espectro elastico 
y la roja al inelastico dividido para el factor B. 





Figura 11 Espectros de analisis para los sismos DBE y MCE 



En base a los perfodos de vibracion se considero T 0 = 2 seg. Valores 
mayores a este perfodo estan relacionados con el sistema de aislacion y menores 
con la estructura. Por lo tanto, para valores de T > T 0 se divide el espectro elastico 
para 5% de amortiguamiento para B. Constantinou etal. (2011). 




Cada aislador tiene un factor de amortiguamiento equivalente. El valor que 
interviene en la ecuacion (8), es el promedio de los factores de amortiguamiento. Se 
encuentra el valor % eq . tanto para el sismo DBE como para el sismo MCE, en base 
a la II y III fase del modelo de Me Vitty y Constantinou (2015). Estos valores ya 
fueron indicados. 

2.9 Nuevo programa para CEINCI-LAB 

En el sistema de computacion CEINCI-LAB, el programa que determina la 
matriz de rigidez en coordenadas de piso, para estructuras con porticos ortogonales 
se denomina matriz_es. Ahora cuando se tienen porticos no ortogonales, como en 
el caso de la estructura analizada se debe utilizar el nuevo programa que han 
elaborado los autores de este artfculo denominado matriz_es1, cuya forma de uso 
se indica a continuation. 



• [KE ]=matriz_es1 (ntot,NP,KL,RT,TETA) 

ntot Numero total de porticos. 

NP Numero de pisos. 
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KL Matriz que contiene las matrices de rigidez lateral de cada portico. 

RT Archivo que contiene todos los vectores r, en todos los pisos. 

TETA Vector que contiene los angulos de cada portico. 

Para el modelo numerico de calculo adoptado, se determino inicialmente la 
determinacion de la matriz de rigidez lateral de los siete porticos. 

Para el ensamblaje de la matriz de rigidez en coordenadas de piso fue 
necesario determinar la distancia (d) perpendicular de los porticos con el Centro de 
Masas CG, descrito en la tabla 8; ademas en el programa que determina la matriz 
de rigidez se ingresa como dato los angulos 0 de cada portico en el mismo orden 
de precedencia en que fueron ingresados los valores de la matriz de rigidez, es decir 
primeramente los porticos en sentido X y posteriormente en sentido Y: 

KT = [KX; KY] = [K44; K47; K48; KY; KZ; Kder; Kiz] 

RT = [RX; RY] = [R1 ; R2; R3; RY; RZ; RD; Rl] 

TETA= [0; 0; 0; 90; 90; 135; 45] 

b = [bx; by; bz]; donde: bx = [1 ; 1; 1]; by = [0]; bz = [0] 

A fin de automatizar el calculo se describe el programa de CEINCI-LAB 
utilizado para el calculo en el sentido X en DBE, para el analisis en el sentido X, Y, 
para sismo DBE y MCE, se sigue el mismo procedimiento: 

H Analisis sismico de bloque 8 en sentido X ante sismo DBE 
% Dr. Roberto Aguiar; Edison Morales; Rodrigo Mejia 
% Diciembre de 2015 

% 

load K44 
load K47 
load K48 
load KY 
load KZ 
load Kder 
load Kiz 

% Matriz de distancias 

% Porticos en sentido Y 

RT* [-6.25 -6.289 -6.53] ;RY*diag(RY) ;% Portico Y 
RZ* [6.25 6.211 5.97] ;RZ*diag(RZ) ;% Portico Z 
%Porticos en sentido X 

Rl* [ 16 . 994 17.883 17 . 723] ; Rl*diag (Rl) ; % Portico 44 

R2- [-1.556 -.667 - . 827] ; R2*diag (R2 ) ; % Portico 47 

R3= [-13 . 906 -13.017 -13 . 177] ; R3*diag (R3 ) ; % Portico 48 

RD* [14.36 13.746 14 . 028] ; RD=diag ( RD) ; % Portico Derecho 

Rl* [-14 . 36 -13.746 -14 . 028] ; RI*diag (Rl) ; % Portico Izquierdo 

nnnnn 

KX* [K44; K47; K48] ; 

KY=[KY;KZ; Kder; Kiz] ; 

KT* [KX; KY] ; %Sent ido X para Y cambiar orden dentro de corchetes 



RX=[R1;R2;R3] ; 

RY* [ RT; RZ ; RD ; Rl] ; 

RT* [RX; RY] ;% Sent ido X para Y cambiar orden dentro de corchetes 
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% [KE] =matr iz_esl (ntot ,NP, KT, RT, TETA) 

% 

% TETA Vector que contiene los angulos de cada Portico 
% ntot Nuinero total de porticos 

% NP Nuinero de Pisos 

% KT= Matriz que contiene las matrices de rigidez lateral de todos los 

% porticos. Priroero los de sentido X. 



RT= 



Vector que contiene la distancia del portico al centro de masa, 
de cada uno de los porticos, con signo, positivo antihorario. Un 
solo valor por portico 

Archivo de datos que contiene vectores r, en todos los pisos. 
Para estructuras regulares se colocara RT=r , pero en estructuras 
irr ^1 valor de R se debe indicar piso por piso en la diagonal. 



TETA= [ 0 ; 0 ; 0 ; 90 ; 90 ; 13 5 ; 45] ; ntot=7;NP=3 ; 
[KE, A] ®matr iz_esl (ntot,NP, KT, RT, TETA) 

% Matriz de Masas 

ml=37.26;m2=56.98;m3=45. 12; 

Jl=4053 . 09; J2=7717 . 68; J3=5940 . 41; 

ME= [ml 00000000 
0 m2 0000000 
m3 0 0 0 0 0 
0 ml 0 0 0 0 
0 m2 0 0 0 
0 0 m3 0 0 
0 0 0 J1 0 



0 0 0 0 J2 



0 0 0 0 0 J3] ; 



[T, f i, OM] =orden eig (KE, ME) ; 



% Factores R y B de modificacion de espectro 

zeda=0.2028; ^Factor de amort iguamiento de los aisladores promedio 
B= ( zeda/0 . 05) *0 . 3 ; H Para roodificar el espectro en rango de aisladores 
R=l; f ip=l . 0; f ie=l . 0; *Para modificar espectro en rango estructura 
to=1.5; %Periodo a partir del cual se reduce espectro por B 
[Ad] =espectro_necl5_mod_DE_l (R, f ip, f ie, B, to, T) 
bx=[l;l;l] ; cer=zeros (NP, 1) ;b= [bx; cer ; cer] ; 

*Factores de participacion modal 
□ for i=l : 3 *NP 

gama ( ; , i) =abs ( ( (f i ( : , i) ) ' *ME*b) / < (f i (:,!))• *ME* (f i ( : , i) ) ) ) ; 

^ end 

yama 

na=3*NP; 

[qt] ®desplazamientos_modales_CQC (T, f i, Ad, gama, na, OM, zeda) %desp elasticos 
qine=R*f ip*f ie^qt; 
qt ine*qine ( 1 : NP ) 

[ Derv] ®der iva ( Y, qt ine ) 

[Vt, Ft] *f uer zas_roodales_CQC (ME, f i, Ad, gama, na, OM, zeda) 

[Fta, Ap] =f uer zas modales aceleracion (ME, f i, Ad, gama, na) 
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2.10 Propiedades dinamicas 

En la tabla 1 0 se indica el periodo y el valor de los modos de vibracion, para 
el sismo DBE y en la tabla 1 1 para el sismo MCE. La estructura tiene 9 grados de 
libertad. 



Tabla 10 Propiedades dinamicas para sismo DBE 



DBE 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


1 


2 


3 


4 


5 


6 


7 


8 


9 


Periodo 


3.78 


3.28 


2.36 


0.47 


0.31 


0.24 


0.06 


0.03 


0.03 




- 0.074 


- 0.003 


0.031 


- 0.134 


0.024 


- 0.043 


0.000 


- 0.000 


0.000 




- 0.076 


- 0.002 


0.040 


0.045 


- 0.008 


0.015 


- 0.004 


0.085 


0.022 




- 0.076 


- 0.002 


0.039 


0.045 


- 0.008 


0.015 


0.004 


- 0.108 


- 0.028 




0.022 


- 0.082 


- 0.001 


0.013 


0.136 


0.040 


0.001 


- 0.000 


0.000 


MODOS DE 
VIBRACION 


0.024 


- 0.086 


- 0.001 


- 0.004 


- 0.047 


- 0.014 


- 0.086 


- 0.008 


0.018 




0.016 


- 0.086 


- 0.002 


- 0.005 


- 0.049 


- 0.012 


0.108 


0.010 


- 0.023 




- 0.003 


- 0.000 


- 0.007 


- 0.004 


- 0.003 


0.013 


- 0.000 


0.000 


- 0.000 




- 0.003 


0.000 


- 0.007 


0.002 


0.001 


- 0.003 


0.002 


- 0.002 


0.007 




- 0.003 


0.000 


- 0.007 


0.002 


0.001 


- 0.004 


- 0.002 


0.002 


- 0.009 



Tabla 11 Propiedades dinamicas para sismo MCE 



MCE 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


1 


2 


3 


4 


5 


6 


7 


8 


9 


Periodo 


3.79 


3.30 


2.38 


0.47 


0.30 


0.24 


0.06 


0.03 


0.03 




- 0.074 


- 0.002 


0.031 


- 0.134 


0.024 


- 0.043 


0.000 


- 0.000 


0.000 




- 0.076 


- 0.002 


0.040 


0.045 


- 0.008 


0.015 


- 0.004 


0.085 


0.022 




- 0.077 


- 0.002 


0.039 


0.045 


- 0.008 


0.015 


0.004 


- 0.108 


- 0.028 




0.002 


- 0.082 


- 0.001 


0.013 


0.136 


0.036 


0.001 


- 0.000 


0.000 


MODOS DE 
VIBRACION 


0.002 


- 0.085 


- 0.000 


- 0.004 


- 0.047 


- 0.014 


- 0.086 


- 0.008 


0.018 




0.001 


- 0.085 


- 0.002 


- 0.005 


- 0.049 


- 0.012 


0.108 


0.010 


- 0.023 




- 0.003 


- 0.000 


- 0.007 


- 0.004 


- 0.003 


0.013 


- 0.000 


0.000 


- 0.000 




- 0.003 


0.000 


- 0.007 


0.002 


0.001 


- 0.004 


0.002 


- 0.002 


0.007 




- 0.003 


0.000 


- 0.007 


0.002 


0.001 


- 0.004 


- 0.002 


0.002 


- 0.009 






Analisis sfsmico espacial y piano del Bloque Estructural 8 



109 



2.1 1 Resultados del analisis sismico en sentido X 

En la figura 10 se indico los porticos en sentido X, corresponden a los que 
tienen menos longitud y se hallan en la direction del voladizo. 

En la tabla 12 se presenta el factor de participacion modal r en los modos 
de vibration y las respectivas aceleraciones espectrales A d para el sismo DBE; en 
la tabla 1 3 se indican los desplazamientos q que se hallan luego de aplicar el criterio 
de combinacion modal CQC y las aceleraciones de piso A p . 

Los valores de la tabla 14 son similares a los de la tabla 12 pero para el 
sismo MCE y los resultados de la tabla 1 5 son similares a los de la tabla 1 3. 



Tabla 12 Factores de participacion modal y aceleraciones espectrales para sismo 

DBE en sentido X 





Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 




1 


2 


3 


4 


5 


6 


7 


8 


9 


r 


10.57 


0.33 


5.23 


0.35 


0.06 


0.10 


0.00 


0.000 


0.000 


A d 

(m/s 2 ) 


1.143 


1.320 


1.832 


11.666 


11.666 


11.666 


8.514 


7.068 


6.656 



Tabla 13 Desplazamie nto y aceleraciones de piso para sismo DBE en sentido X 





Nivel + 0.04 


Nivel + 4.94 


Nivel + 9.88 


q (cm) 


31.02 


31.64 


31.90 


F (T) 


16.72 


20.75 


16.18 


A p (m/s 2 ) 


1.09 


1.02 


1.02 



Tabla 14 Factores de participacion modal y aceleraciones espectrales para sismo 

MCE en sentido X 





Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 




1 


2 


3 


4 


5 


6 


7 


8 


9 


r 


10.58 


0.30 


5.22 


0.35 


0.06 


0.10 


0.00 


0.00 


0.00 


A d 

(m/s 2 ) 


1.632 


1.873 


2.592 


16.332 


16.332 


16.332 


11.919 


9.895 


9.319 



Tabla 15 Desplazamie nto y aceleraciones de piso para sismo MCE en sentido X. 





Nivel + 0.04 


Nivel + 4.94 


Nivel + 9.88 


q (cm) 


44.48 


45.42 


45.80 


F (T) 


23.72 


29.72 


23.11 


A p (m/s 2 ) 


1.55 


1.45 


1.45 
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DBE MCE 



m3 




O) 




Figura 12 Resultado de desplazamientos en sentido X, para sismo DBE y MCE 



2.12 Resultados del analisis sismico en sentido Y 



Para el ensamblaje de la matriz de rigidez en coordenadas de piso, fueron 
ingresados los valores en el orden de precedencia en que fueron ingresados los 
valores de la matriz de rigidez, es decir primeramente los porticos en sentido Y y 
posteriormente en sentido X: 

KT = [KY; KX] = [KY; KZ; Kder; Kiz; K44; K47; K48] 

RT = [RY; RX] = [RY; RZ; RD; Rl; R1 ; R2; R3] 

TETA= [90; 90; 135; 45; 0; 0; 0] 

b = [bx; by; bz]; donde: bx = [0]; by = [1 ; 1 ; 1]; bz = [0] 



Cabe mencionar que la matriz rigidez en coordenadas de piso, tanto en 
sentido X como en sentido Y es la misma. 



Los resultados para el sismo DBE se indican en las tablas 16 y 17; para el 
sismo MCE en las tablas 1 8 y 1 9. 



Tabla 16 Factores de participation modal y aceleraciones espectrales para sismo 

DBE en sentido Y 





Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 




1 


2 


3 


4 


5 


6 


7 


8 


9 


r 


0.43 


11.79 


0.20 


0.04 


0.17 


0.03 


0.00 


0.000 


0.000 


A d 

(m/s 2 ) 


1.379 


1.530 


2.027 


1 1 .666 


11.666 


11.666 


8.514 


7.068 


6.656 



Tabla 17 Desplazamie nto y aceleraciones de piso para sismo DBE en sentido Y 





Nivel + 0.04 


Nivel + 4.94 


Nivel + 9.88 


q (cm) 


44.01 


45.54 


45.56 


F (T) 


53.96 


83.96 


58.77 
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A p (m/s 2 ) 


1.51 


1.54 


1.54 



Tabla 18 Factores de participation modal y aceleraciones espectrales para sismo 

MCE en sentido Y 





Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 


Modo 




1 


2 


3 


4 


5 


6 


7 


8 


9 


r 


0.38 


11.80 


0.17 


0.04 


0.17 


0.03 


0.00 


0.00 


0.00 


A d 

(m/s 2 ) 


1.630 


1.800 


2.383 


16.332 


16.332 


16.332 


11.919 


9.895 


9.319 



Tabla 19 Desplazamie nto y aceleraciones de piso para sismo MCE en sentido Y. 





Nivel + 0.04 


Nivel + 4.94 


Nivel + 9.88 


q (cm) 


52.26 


54.06 


54.10 


F (T) 


64.21 


101.98 


77.56 


A p (m/s 2 ) 


1.79 


1.82 


1.82 



Figura 13 Resultado 




3. ANALISIS SISMICO PLANO 

El Bloque Estructural 8 tiene 9 aisladores, que han sido identificados en la 
figura 14, donde ademas se detalla las areas cooperantes que aportan a los 
aisladores de los porticos 44 y 48 en cada nivel. 
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Nivel + 0.04 



Nivel + 9.88 

Figura 14 Areas cooperantes que cagan a los aisladores 

A manera de ejemplo se indica en las tablas 20 y 21 el calculo de las cargas 
que Megan a los aisladores 44 V, 48 Y. 
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Tabla 20 Calculo de la Carga que Mega al aislador 44 V 



Aislador 44 V 


Nivel 


Area Cooperante 


Carga Vertical 


W 


(m 2 ) 


(T/m 2 ) 


(T) 


+ 0.04 


31.34 


0.86 


27.03 


+ 4.94 


44.70 


0.86 


38.55 


+ 9.88 


44.70 


0.66 


29.61 




Total 


95.19 



Tabla 21 Calculo de la Carga que Mega al aislador 48 Y 



Aislador 48 Y 


Nivel 


Area Cooperante 


Carga Vertical 


W 


(m 2 ) 


(T/m 2 ) 


(T) 


+ 0.04 


48.61 


0.86 


41.92 


+ 4.94 


80.07 


0.86 


69.06 


+ 9.88 


80.07 


0.66 


53.04 




Total 


164.03 



3.1 Rigidez efectiva de los aisladores 

En el analisis sismico espacial, se presento el aislador FPT8833/1 2-1 2/8-6, 
que se colocara en el Bloque Estructural 8, el mismo que en su geometrfa y 
caracteristicas de diseno rigen tambien para el analisis sismico piano. A manera de 
ejemplo se indica en las tablas 22 y 23 el calculo de las rigideces efectivas de los 
aisladores en los Porticos 44 y 48. 



Tabla 22 Calculo de la rigidez efectiva de los aisladores Portico 44 Bloque 8 



Tabla resumen del Portico 44 Bloque 8 


Aislador 


w 

(T) 


Fmax 2 
(T) 


q2 

(m) 


q3 

(m) 


ss'e 

CD ^ 

SC t 


Kef3 

(T/m) 


44Y 


95.19 


18.022 


0.537 


0.566 


33.586 


31.824 


44Z 


109.31 


20.696 


38.568 


36.545 



Tabla 23 Calculo de la rigidez efectiva de los aisladores Portico 48 Bloque 8 



Tabla resumen del Portico 48 Bloque 8 


Aislador 


w 

(T) 


Fmax 2 
(T) 


q2 

(m) 


q3 

(m) 


CD ^ 

sc H 


Kef3 

(T/m) 


48Z 


164.03 


31.056 


0.537 


0.566 


57.876 


54.840 


48Y 


164.03 


31.056 


57.875 


54.840 
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3.2 Geometria de los Porticos 

En la figura 15 y 16 se muestra la geometria de todos los porticos 44 y 48 
para el analisis sismico piano. 



2.070 tf/m 1 880 tf/m 



0.880 tf/m 

sTTT^TTTi 


innuinuj.uunuujuunu.uiuj 


lUJUUJ J J 




2.700 tf/m 


2.450 tf/m 


1.150 tf/m 


.JUJJUJJJUUJJJUJJJUJUJJJUUU. 


.UJUJJJJJ 




2.700 tf/m 






.UU1UU1UU1UUU1U1UU1UUU. 




Aislador 

77 


77 77 . 


Aislador 

77 



|-£— 2 65 m > j C 12 50 m 3.61 m 

Figura 15 Portico 44 analisis sismico piano 



3.000 tf/m 



1 540 tf/m 

^TTTTTTTTi 


UJJJJJJJJJJJJJJJJJJJUJJUJJ UJUJJJJJ 


1 540 tf/m 

rTTTTTTTTT. 


2.010 tf/m 

xtttttttt: 


3 900 tf/m 


2.010 tf/m 

TTTTTTTTl 


.jjjjjjjnjjinunijjjjjjnjjjjjujjT 


Aislador 

77 . 


3 900 tf/m 


Aislador 

77 


.JJJJUJUJJJJJJJJJJJJJJJJUJUJJUJJJ. 


77 77 



3 10 m 12.53 m 3.10 m 

Figura 16 Portico 48 analisis sismico piano 



Para determinar la carga distribuida equivalentes en los respectivos niveles 
de cada portico se realizo mediante el calculo de mosaico de cargas por nivel en 
cada portico como se muestra en la figura 1 7. 
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Figuras 17 Mosaico de Cargas para los porticos 44 y 48 
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En las tablas 24 y 25, se determinan el calculo de cargas distribuidas y 
masas los porticos 44 y 48 para el analisis sfsmico piano. 

Tabla 24 Cargas distribuidas y masas del portico 44 



Portico 44 


Nivel 


Carga 

(T/m 2 ) 


Vano 1 
(m) 


Area 

1 

(m 2 ) 


CD 1 
(T/m) 


Masa 1 
(T*s 2 /m) 


Vano 2 
(m) 


Area 2 
(m 2 ) 


CD 2 
(T/m) 


Masa 2 
(T*s 2 /m) 


Vano 

3 

(m) 


Area 3 
(m 2 ) 


CD 3 
(T/m) 


Masa 3 
(T*s 2 /m) 


Masa 

Total 

(TV/m) 


9.88 


0.66 


2.65 


3.53 


0.88 


0.24 


12.50 


39.07 


2.07 


2.64 


3.61 


10.25 


1.88 


0.69 


3.57 


4.94 


0.86 


2.65 


3.53 


1.15 


0.31 


12.50 


39.07 


2.70 


3.43 


3.61 


10.25 


2.45 


0.90 


4.65 


0.04 


0.86 


0 


0.00 


0.00 


0.00 


12.50 


39.07 


2.70 


3.43 


3.61 


0.00 


0.00 


0.00 


3.43 



Tabla 25 Cargas distribuidas y masas del portico 48 



Portico 48 


Nivel 


Carga 

(T/m 2 ) 


Vano 1 
(m) 


Area 

1 

(m 2 ) 


CD 1 
(T/m) 


Masa 1 
(T*s 2 /m) 


Vano 2 
(m) 


Area 2 
(m 2 ) 


CD 2 
(T/m) 


Masa 2 
(T*s 2 /m) 


Vano 

3 

(m) 


Area 3 
(m 2 ) 


CD 3 
(T/m) 


Masa 3 
(T*s 2 /m) 


Masa 

Total 

(T*s 2 /m) 


9.88 


0.66 


3.10 


7.22 


1.54 


0.49 


12.26 


55.45 


3.00 


3.74 


3.10 


7.22 


1.54 


0.49 


4.72 


4.94 


0.86 


3.10 


7.22 


2.01 


0.64 


12.26 


55.45 


3.90 


4.88 


3.10 


7.22 


2.01 


0.64 


6.15 


0.04 


0.86 


0 


0.00 


0.00 


0.00 


12.26 


55.45 


3.90 


4.88 


0 


0.00 


0.00 


0.00 


4.88 



3.3 Propiedades dinamicas Portico 44 

En la tabla 26 se indica el periodo y el valor de los modos de vibracion, para 
el sismo DBE y en la tabla 27 para el sismo MCE. La estructura tiene 3 gdl. 

Tabla 26 Propiedades dinamicas para sismo DBE 



DBE 


Modo 1 


Modo 2 


Modo 3 


Periodo 


2.57 


0.270 


0.017 


MODOS DE 
VIBRACION 


0.284 


-0.459 


0.001 


0.297 


0.184 


0.305 


0.297 


0.182 


-0.398 



Tabla 27 Propiedades dinamicas para sismo MCE 



MCE 


Modo 1 


Modo 2 


Modo 3 


Periodo 


2.637 


0.270 


0.017 


MODOS DE 
VIBRACION 


0.284 


-0.459 


0.001 


0.297 


0.184 


0.305 


0.297 


0.183 


-0.399 
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Tabla 28 Factores de participation modal y aceleraciones espectrales para sismo 





DBE en sentido X 




Modo 


Modo 


Modo 




1 


2 


3 


r 


3.413 


0.068 


0.000 


A d 

(m/s 2 ) 


1.68 


11.666 


5.845 



Tabla 29 Desplazamie nto y aceleraciones de piso para sismo DBE en sentido X 





Nivel + 0.04 


Nivel + 4.94 


Nivel + 9.88 


q (cm) 


27.29 


28.53 


28.52 


F ( T ) 


6.01 


8.12 


6.23 


A p (m/s 2 ) 


1.67 


1.71 


1.71 



Tabla 30 Factores de participation modal y aceleraciones espectrales para sismo 





MCE en sentido X 




Modo 

1 


Modo 

2 


Modo 

3 


r 


3.41 


0.06 


0.00 


A d 

(m/s 2 ) 


2.34 


16.33 


8.18 



Tabla 31 Desplazamie nto y aceleraciones de piso para sismo MCE en sentido X. 





Nivel + 0.04 


Nivel + 4.94 


Nivel + 9.88 


q (cm) 


40.06 


41.77 


41.77 


F ( T ) 


8.34 


11.28 


8.66 


A p (m/s 2 ) 


2.32 


2.38 


2.34 



DBE MCE 




Figura 18 Resultado de desplazamientos en sentido X, para sismo DBE y MCE 



118 



Roberto Aguiar, Rodrigo Mejia y Edison Morales 



3.4 Propiedades dinamicas Portico 48 

En la tabla 31 se indica el periodo y el valor de los modos de vibracion, para 
el sismo DBE, en la tabla 32 para el sismo MCE. 

Tabla 32 Propiedades dinamicas para sismo DBE 



DBE 


Modo 1 


Modo 2 


Modo 3 


Periodo 


2.382 


0.316 


0.020 


MODOS DE 
VIBRACION 


0.240 


-0.384 


0.001 


0.257 


0.161 


0.265 


0.257 


0.160 


-0.347 



Tabla 33 Propiedades dinamicas para sismo MCE 



MCE 


Modo 1 


Modo 2 


Modo 3 


Periodo 


2.44 


0.316 


0.020 


MODOS DE 
VIBRACION 


0.241 


-0.383 


0.001 


0.257 


0.162 


0.265 


0.257 


0.160 


-0.347 



Tabla 34 Factores de participacion modal y aceleraciones espectrales para sismo 





DBE en sentido X 




Modo 


Modo 


Modo 




1 


2 


3 


r 


3.97 


0.13 


0.00 


A d 

(m/s 2 ) 


1.82 


11.66 


6.02 



Tabla 35 Desplazamie nto y aceleraciones de piso para sismo DBE en sentido X 





Nivel + 0.04 


Nivel + 4.94 


Nivel + 9.88 


q (cm) 


24.86 


26.64 


26.64 


F (T) 


9.45 


11.83 


9.08 


A p (m/s 2 ) 


1.82 


1.87 


1.87 



Tabla 36 Factores de participacion modal y aceleraciones espectrales para sismo 

MCE en sentido X 





Modo 


Modo 


Modo 




1 


2 


3 


r 


3.97 


0.12 


0.00 


A d 

(m/s 2 ) 


2.53 


16.33 


8.43 
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Tabla 37 Desplazamie nto y aceleraciones de piso para sismo MCE en sentido X. 





Nivel + 0.04 


Nivel + 4.94 


Nivel + 9.88 


q (cm) 


36.50 


38.98 


38.97 


F ( T ) 


15.67 


19.65 


15.09 


A p (m/s 2 ) 


2.53 


2.60 


2.60 



DBE MCE 




Figura 19 Resultado de desplazamientos en sentido X, para sismo DBE y MCE 

4. COMENTARIOS Y CONCLUSIONES 

Se ha realizado un analisis sismico espacial y piano del Bloque Estructural 
8 de la UFA-ESPE que tiene porticos inclinados y una geometria especial en planta 
con voladizos en los pisos superiores. Se ha presentado con bastante detalle el 
calculo realizado en los porticos inclinados y se ha incorporado un nuevo programa 
al sistema de computacion CEINCI-LAB que sirve para encontrar la matriz de rigidez 
en coordenadas de piso en estructuras cuyos porticos no son ortogonales. 

En forma cuantitativa se puede decir que los resultados hallados en el 
analisis sismico espacial son comparables con los del analisis sismico piano, en el 
sentido de que muestran que es el aislador el que mas se desplaza y el 
desplazamiento relativo de los porticos con respecto al sistema de aislacion es muy 
bajo; eso se aprecio con los dos modelos. En forma cualitativa hay diferencia en los 
resultados y la causa principal es que en el analisis sismico piano no se considera 
la torsion, en cambio en el analisis sismico espacial si se tiene en cuenta la torsion. 
Otra de las causas es la determinacion de las masas. 

Se aspira haber aportado al desarrollo de la Ingenieria Sismica, al presentar 
en forma bastante detallada el analisis sismico de una estructura compleja con 
aisladores sismicos, que dicho sea de paso que a pesar de su complejidad la 
superestructura experimenta desplazamientos laterales relativos al desplazamiento 
del aislador, muy bajos, lo que demuestra la bondad de tener estructuras sobre 
aisladores sismicos. 
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